As part of the effort to identify and classify the ob- 
jects responsible for this X-ray emission, we have ob- 
tained optical spectra of the brightest optical counterparts 
of the identified X-ray sources, using the 1.3m Skinakas 
Telescope. Most of these objects are foreground star can- 
didates. The purpose of the present study is to confirm this 
identification and to explore the compatibility between the 
optical spectral classification and the observed X-ray prop- 
erties of the sources. We have obtained optical spectra for 
the 14 brightest optical counterparts of X-ray sources iden- 
tified by XMM-Newton in the direction of M 31 and for 
21 optical counterparts in the direction of M 33, using the 
1.3m Skinakas telescope in Crete, Greece. All of the M 31 
\^ ■ sources and all but one of the M 33 sources were confirmed 
[ to be foreground stars, of spectral types between A and 
■ M.One of the stars is a late M dwarf with emission, 
a flare star, also displaying strong X-ray variability. One 
of the M 33 sources (lying within the D25 ellipse) corre- 
sponds to a previously known background galaxy, LEDA 
5899. 
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Abstract. Recent surveys of the Local Group spiral Galaxies M 31 and M 33 with XMM-Newton yielded a large 
number of X-ray sources. 



1. Introduction 

Recent surveys of the Local Group spiral Galaxies M 31 
(Pietsch et al. 2005a, hereafter PFH2005) and M 33 
(Pietsch et al. 2004, hereafter PMH2004) with XMM- 
Newton yielded a large number of X-ray sources. With 
their moderate Galactic foreground absorption, both 
galaxies are ideal for studying the X-ray source popula- 
tion and diffuse emission in nearby spiral galaxies. 

The Andromeda galaxy (M 31) is a massive SA(s)b 
galaxy, located at a distance of 780 kpc (Holland 1998) 
and seen under an inclination of 78° (de Vaucouleurs et 
al. 1991). 

As described in detail in PFH2005, a total of 856 X- 
ray sources were detected in an area of 1.24 square de- 
grees. Within M 31, 21 supernova remnants (SNR) and 
23 SNR candidates were detected, as well as 18 super-soft 
source candidates, 7 X-ray binaries and 9 X-ray binary 
candidates, 27 globular cluster sources and 10 globular 
cluster source candidates, which most likely are low mass 
X-ray binaries within the cluster. From 567 hard sources, 
some are expected to be X-ray binaries or Crablike SNRs 
within M 31, and the rest, background AGN. As well as 
the sources within M 31, there were also foreground and 
background sources: six sources were identified as fore- 
ground stars and 90 as foreground star candidates, one 
as a BL Lac type active galactic nucleus (AGN) and 36 
as AGN candidates. One source coincided with the Local 
Group galaxy M 32, one with a background galaxy cluster 
and another is a background galaxy cluster candidate. 

M 33, an Sc galaxy, is located at a distance of 795 
kpc (van den Bergh 1991) and is seen under a relatively 
low inclination of 56° (Zaritsky et al. 1989). PMH2004 de- 
tected a total of 408 sources in a 0.8 square degree field 
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combining the counts of all EPIC instruments. About half 
of the sources lie within M 33; the rest are background 
AGNs or foreground stars. Cross-correlation with archival 
data led to the identification of 5 foreground stars and 
30 foreground star candidates, as well as 12 AGN candi- 
dates. Within M 33, 21 supernova remnants (SNR) and 
23 SNR candidates were identified, as well as 5 super-soft 
sources and 2 X-ray binaries. From the remaining 267 X- 
ray sources classified as "hard" , those within M 33 may be 
low or high mass X-ray binaries, or Crab-like SNRs, while 
a significant fraction of them are expected to be back- 
ground AGN. Misanovic et al. (2005) analyzed the indi- 
vidual observations of the XMM-Newton data of M 33, de- 
tecting 39 new sources and improving the positions of an- 
other 311. They also studied the variability of the sources 
on time scales of hours to months or years, their spectral 
characteristics and classification. For example, using the 
detected variability, they were able to classiiy 8 new X-ray 
binary candidates. 

Pietsch et al. (2005b) searched for X-ray counterparts 
of optical novae detected in M 31 and M 33 and discovered 
21 X-ray counterparts for novae in M 31 - mostly identified 
as supersoft sources by their hardness ratios - and two in 
M 33. 

As part of the effort of identifying and classifying the 
objects responsible for the observed X-ray emission in the 
direction of M 31 and M 33, we have obtained optical spec- 
tra of the brightest optical counterparts of the identified 
X-ray sources, using the 1.3m Skinakas Telescope. These 
brightest optical counterparts were taken from PFH2005 
for M 31 and from PMH2004 for M 33. Most of these ob- 
jects are foreground star candidates, as also mentioned in 
the corresponding reference papers. The purpose of the 
present study is to confirm this and to explore the com- 
patibility between the optical spectral classification and 
the observed X-ray properties of the sources. One of the 
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bright optical counterparts that we have observed is the 
known galaxy LED A 5899, which was also part of the CFA 
redshift survey (Huchra et al. 1999). 

In the next section, we describe the optical observa- 
tions; in Section 3 we describe the data reduction. In 
Section 4. the method followed for the classification of the 
spectra obtained is presented and in Section 5 we discuss 
the results, including individual objects of interest. 

2. Optical observations 

The optical observations used in this study were carried 
out during three observing runs, on September 6, 2003, on 
October 18-22, 2003 and on October 10-11, 2004, using the 
1.3-m Ritchcy-Crcticn telescope at Skinakas Observatory, 
located on the island of Crete (Greece). 

The telescope was equipped with a 2000x800 ISA SITe 
CCD camera and a 1302 lines/mm grating, giving a nom- 
inal dispersion of 1.04 A/pixel, and a wavelength coverage 
from 4748 A to 6828 A. The spectral region was selected 
so as to include both Hq and Hp spectral lines. For the 
three spectra that were obtained in September 2003, the 
wavelength region was somewhat different, shifted further 
to the blue (from 3900 A to 6000 A). One of these three 
spectra (of source number 3 in M 33) was also observed 
in October 2003, with the wavelength coverage used for 
the bulk of the spectra. As will be discussed later, the two 
spectra obtained for this object served as a test of the 
compatibility between the spectral classes obtained in the 
two different wavelength regions. 

Exposure times ranged from 600 to 7200 s, depend- 
ing on the object's magnitude and on seeing and weather 
conditions. In most cases, two or more exposures were ob- 
tained per object. Each object observation was followed 
by an arc calibration exposure (CuAr). 

A total of 18 standards, of spectral types ranging from 
09.5 to M6, were also observed with the same configura- 
tion as the targets. 

In the first observing run three objects in the direc- 
tion of M 33 were observed (6/9/2003). In the second run 
another 19 targets in the direction of M 33 were observed 
(one of them repeated from the first run, so a total of 21 
targets in the direction of M 33 were observed), while in 
the third, spectra for 14 objects in the direction of M 31 
were obtained. The first 6 columns of Table 1 list the de- 
tails of the observations. 

1. Column 1 gives the X-ray source number of the tar- 
get as it appears in the corresponding X-ray catalogue 
paper, i.e., in PFH2005 for M 31 and in PMH2004 for 
M 33. There is one exception, object 3b in the direc- 
tion of M 33, which was not detected by the automatic 
procedures due to its far off axis location and thus it 
was not included in PMH2004 (see Section 5.2 for more 
details). 

For object 122 in the direction of M 33, the bright 
USNO-Bl star 1210-0020531 given as possible optical 
counterpart in PMH2004 is actually resolved into two 



Fig. 1. X-ray positions (3-sigma error circles in black) of 
PMH2004 sources overlaid on the optical image of Massey 
et al. (2002). The bright stellar counterparts with spectro- 
scopic follow-up of [PMH2004] 297 and 281 are clearly visi- 
ble. The smaller white circles show the improved positions 
and errors from the analysis of individual XMM-Newton 
observations (Misanovic et al. 2005). Note that object 281 
was not detected in this latter study. The figure also shows 
one additional X-ray source in the FoV that was classified 
as 'hard' ([PMH2004] 304). 

stars on deeper 4m Mosaic images of M 33, separated 
by ~6" in the NNE/SSW direction (see Section 5.3 
for more details). Spectra for both of these stars were 
obtained. 

If available, Chandra source names are given in the 

footnotes of the Table. For M 31 only one source 
([PFII2005] 31) has a Chandra identification, while for 
M 33 two sources ([PMH2004] 196 and 200) have also 
been identified in existing Chandra catalogues. 

2. Columns 2& 3 give the co-ordinates of the objects 
(epoch 2000). 

3. Column 4 provides the USNO-Bl identification derived 
from the cross-correlation between the X-ray catalogue 
and the USNO-Bl catalogue. 

4. Columns 5 & 6 list the B and R magnitudes of the 
corresponding USNO object (B2 and R2 in the USNO- 
Bl catalogue). The limiting magnitude of the sam- 
ple of stars observed in M 33 was R=14.75, while in 
M 31 it was 12.43. In neither case is the sample of ob- 
jects observed complete to the corresponding limiting 
magnitude. The brighter limiting magnitude for the 
M 31 stars was dictated by observing time limitations 
only (as three of the nights allocated were lost to poor 
weather conditions). 

5. Column 7 gives the 2MASS identification derived from 
the cross-correlation between the X-ray catalogue and 
the 2MASS catalogue. 

6. Column 8 provides the corresponding J magnitude of 
the 2MASS counterpart (as given in the 2MASS cata- 
logue). 

7. Column 9 gives the date of the Skinakas observation. 

8. Column 10 lists the total exposure time for the specific 
spectrum, while the number of exposures taken for that 
object is provided in the parenthesis following the total 
exposure time. 

Figure 1 illustrates how well the X-ray positions cor- 
relate with the coordinates of the optical counterparts ob- 
served. 

3. Data reduction 

The data reduction was performed using the STARLINK 
Figaro package (Shortridge et al. 2001). The frames were 
bias subtracted, flat flelded and corrected for cosmic ray 
events. The 2-D spectra were subsequently sky subtracted 
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Table 1. Log of objects observed in the direction of M 31 and M 33 and derived spectral types. 
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* X-ray source numbers from PFH2005 (M 31) and PMH2004 (M 33) 
^ Chandra source sl-74, from Williams et al. 2004 

The adopted spectral type for this object is the average of the spectral types obtained from the two observing runs. 
^ This object is a newly discovered source and it was not included in the PMH2004 catalogue. See section 5.2 for details. * far 
off-axis in Chandra observations 1730, not in source list of Grimm et al. (2005) 

^ As explained in Section 5.3, the bright USNO-Bl star 1210-0020531 is resolved into two stars on deeper 4m Mosaic images of 
M 33, separated by ^6" in NNE/SSW direction. The coordinates of these two stars and the corresponding spectral types we 
derived for them, are given in the two lines following the entry for source 122. 

® Chandra source CXO J013341. 8+303848, from Grimm et al. 2005; no USNO-Bl position, 2MASS position given. 
Chandra source CXO J013343.4+304630, from Grimm et al. 2005 
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using the POLYSKY commaiid. A spatial profile was then 
determined for each 2-D spectrum, and the object spectra 
were optimally extracted using the algorithm of Horne 
(1986), with the OPTEXTRACT routine. Arc spectra 
were then extracted from the arc exposures, using exactly 
the same profiles as for the corresponding object spectra. 
The arc spectra were subsequently used to calibrate the 
object spectra. After calibration, spectra from individual 
exposures of the same object were co-added to yield the 
final spectrum. 

This procedure was followed for almost all of the ob- 
servations, with one exception, namely object USNO-Bl 
1210-0020531 in the M 33 list. In this second ob- 

ject was located very close to the target (see section 5.2). 
The resulting spectra could not be extracted separately 
using the optimal extraction procedure described above. 
Due to the slight curvature of the spectra, single linear 
extraction along the lines of the CCD was not possible 
either. Thus we have used a step extraction with suffi- 
cient wavelength overlap between subsequent extractions 
to ensure compatibility of the fluxes. The calibrating arc 
spectra were created in exactly the same way. 

The signal-to-noise ratio (S/N) of the final spectra var- 
ied from 20 to about 350, with an average of 150. 

The same reduction and calibration procedures were 
followed for the 18 spectroscopic standard stars. Most of 
the standard spectra have a S /N ranging between 200 and 
350, while two have somewhat lower S/N (around 130). 

Figure 2 shows six examples of target spectra (black 
lines), all of foreground stars of different spectral types, 
that have been flux-normalized for presentation purposes. 
For comparison, we show in grey the corresponding stan- 
dard star spectrum (of the same spectral type) shifted in 
y by an arbitrary amount. 

4. Spectral classification 

Classification of the obtained spectra was achieved via 
cross-correlation with the 18 standard star spectra. 

Each object spectrum was cross-correlated with each 
standard spectrum, and the height of the correspond- 
ing cross-correlation peak (hereafter, ccp) was recorded. 
These ccp heights were plotted as a function of the spec- 
tral type of the standard, with the maximum of the curve 
yielding the adopted spectral type for the object spec- 
trum. Figure 3 shows a typical example of such a plot, for 
object USNO-Bl 1201-0020193. The highest ccp is about 
0.9, corresponding to an F5 spectral type. For the entire 
sample, the highest ccp, on which the spectral classifica- 
tion was based, ranged from 0.7 to 0.95, with most being 
higher than 0.8. 

In some cases there were two ccp (usually for adjoining 
spectral types in the standard star grid) with comparable 
heights. For these, we adopted the average spectral type 
as best representing the spectral class of the object. 

As an additional check, we inspected the spectra vi- 
sually and compared them against those of the standard 
stars. In very few cases did the spectral type yielded by 



Fig. 3. Results of cross-correlation between the spectrum 
of object USNO-Bl 1201-0020193 with each standard star 
spectrum. On the x-axis, the spectral type of the standard 
star is indicated, while on the y-axis the height of the 
corresponding cross-correlation peak is given. 

Fig. 4. Derived versus standard spectral type for the stan- 
dard stars 



the cross-correlation method described above have to be 
modified, after visual inspection. These modifications were 
always within our estimated error of 0.3 of a spectral type 
(see below). 

The adopted spectral type for each objects is reported 
in the last column of Table 1. 

The accuracy of the spectral classification achieved in 
this manner depends on the fineness of the grid of standard 
spectra used and on the signal-to-noisc ratio of the cross- 
correlated spectra. We estimated the accuracy of the spec- 
tral classification in the following manner: We treated each 
standard star as an object spectrum. We cross-correlated 
it with all of the other standards, and determined its spec- 
tral type using the cross-correlation peak curves, as was 
done for the object spectra. In Fig. 4 we plot the derived 
spectral type for the standards as a function of the refer- 
ence spectral type. A least-squares linear fit yields a slope 
of 0.98 ± 0.06 and a scatter of 0.3 in spectral type. 

The effect of the S/N on the derived spectral type was 
estimated by utilizing the individual exposures of vari- 
ous targets. Following the same procedure as for the co- 
added spectra we derived spectral types and compared 
them to the spectral types adopted for the co-added spec- 
tra. Although the cross-correlation peaks were lower for 
the lower S/N spectra, the maximum value of the cross- 
correlation peaks was equally well defined and no differ- 
ences in spectral classification were noted down to a S/N 
of about 30. There are two spectra for which the S/N 
is significantly lower and for these the spectral classifica- 
tion is less accurate (spectral type followed by the symbol 
"::", in Table 1). Therefore, for almost all of our spectra 
the accuracy of the spectral classification depends on the 
fineness of the standard-star grid and is around 0.3 of a 
spectral type. 

For the three spectra obtained in September 2003, the 
cross-correlation was performed over the common wave- 
length range between these spectra and the standards 
(that were obtained in the subsequent observing runs), 
i.e. over lOOOA. Despite the shorter wavelength range, the 
cross-correlation peaks were very high, similar to those 
obtained for the rest of the spectra. As mentioned earlier, 
one of the stars ([PMH2004] 3 in M 33) was observed in 
both wavelength ranges. The spectral classes derived were 
in excellent agreement (within 0.1 of a spectral type). 
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Fig. 2. Six examples of spectra (flux normalized) of optical counterparts of X-ray sources in the direction of M 31 and 
M 33 (black lines). For comparison, we show in grey the corresponding standard star spectrum (of the same spectral 
type) shifted in y by an arbitrary. The identified lines are H„ (6563A), H/j (4861 A), Fel (5270, 5328, 6495A) and Nal 
(5890A). 



5. Discussion 

All the objects observed, except for [PMH2004] 358 in the 
direction of M 33 which is a known galaxy (LEDA 5899) 
and which is discussed separately below, arc consistent 
with being foreground stars, the spectral types of which 
range from A5 to M5. Figure 5 shows the frequency dis- 
tribution of the derived spectral typos. 

In Table 2, we list some of the X-ray properties of the 
sources (as given in PMH2004 and PFH2005), namely the 
X-ray flux, fx , the logarithm of the ratio of the X-ray flux 
to the optical flux, log{fx/ fopt) and the hardness ratios 
HRl and HR2. 

The fluxes fx were calculated from the 0.2-4.5 keV 
XID band. The ratios log{fxl fopt) were calculated using 
the formula log{fx/ fopt) = log{fx) + (B + R)/{2 x 2.5) + 
5.37(Maccacaro et al. 1988) and the USNO-bl B and R 
magnitudes given in Table 1. The hardness ratios were 
derived from the count rates in the energy bands R\ = 
0.2 - 0.5 keV, R2 = 0.5 - 1.0 keV, and i?3 = 1.0 - 2.0 keV, 
according to the formulae HRl = {R2 — Ri)/{R2 + Ri) 
and HR2 = {R3 - R2)/{Rz + R2), as in PMH2004. 

The X-ray fluxes span two orders of magnitude, rang- 
ing from 1.5e-15 mW/m^ to 1.47e-13 mW/m^ . The values 
of log{fx/ fopt) range from -5.9 to -2.5, with an average of 
-4.2 ±0.8. HRl takes values between 0.16 and 0.94, with 
an average of 0.5 ± 0.2, while HR2 ranges from -0.93 to 
0.32, with an average of —0.5 ± 0.3. There is no clear cor- 
relation between any of the X-ray properties listed here 
and the derived spectral types of the optical counterparts. 
However, the object with the lowest X-ray to optical flux 
ratio is the only early type star in the sample (of A type). 

All of the objects included in the present study, ex- 
cept for the known galaxy LEDA 5899, satisfy the basic 
criterion for being foregroimd stars, i.e. log{fx/ fopt) < —1 
(Maccacaro et al. 1988 ). They also satisfy the additional 
criterion described in PFH2005, HR2-EHR2<0.3 (where 
EHR2 is the error in the value of HR2, also given in Table 
2). In Figure 6 we plot the hardness ratio HR2 against the 
X-ray to optical flux ratio [log{fx/ fopt)], marking the dif- 
ferent spectral classes in different colours. It is clear that 
all spectral classes in our sample (with enough stars to 
make the comparison meaningful) occupy the same locus 
on the log{fx/ fopt)-tl^2 plane. However, there are three 
objects, i.e. [PFH2005] 464 in the direction of M 31 and 
[PMH2004] 206 and [PMH2004] 337 in the direction of 
M 33, that have hardnc;ss ratios too high compared to the 
rest of the sample (although they do satisfy the previously 
mentioned criterion for foreground stars). Further inspec- 
tion of these sources showed that two of them may have 
overestimated hardness ratios due to their proximity to 
bright hard sources (see discussion of individual objects 



Fig. 5. Frequency distribution of spectral types of the ob- 
served bright stellar counterparts of X-ray sources in M 31 
and M 33 

Fig. 6. X-ray properties of the observed sample. Different 
spectral types are marked in different colours. 

Fig. 7. X-ray contours from the (0.2 1.0) keV combined 
image of PMH2004 overlaid on the optical R image of 
Massey et al. The contours are at levels 3, 4 and 6 x 10~^ 
cts pixel""*^, and indicate the position of additional X- 
ray emission coinciding with the USNO-Bl source 1204- 
0019500. The position on a pn CCD gap prevented the 
detection of this emission as an X-ray source. The nearby 
source [PMH2004] 3 is shown with a box whose size cor- 
responds to the 3 sigma error. 

below). Further discussion of the three sources follows in 
paragraphs 5.1, 5.5 and 5.7. 

5.1. [PFH2005] 464, in the M 31 field (USNO-Bl 
1311-0013521) 

This object has a hardness ratio (HR2) value that is high 
compared to other stars in the sample with similar X- 
ray to optical flux ratios (as seen in Fig. 6). The optical 
spectrum of the star shows no indication of emission lines 
that might indicate flaring activity. However, the X-ray 
source itself is quite faint and relatively close to two bright 
hard sources [PFH2005] 457 (fx = 3.9e - 13 ± 5e - 15, 
HR2 = 0.46 ± 0.01) and [PFH2005] 463 {fx = 2.7e - 13 ± 
4e — 15, HR2 = 0.81 ±0.02). Both sources are at a distance 
of 1.52 arcmin from [PFH2005] 464 and have higher X-ray 
fluxes by two orders of magnitude. This may disturb the 
X-ray detection program in the determination of the hard 
count rate and lead to a harder HR2. 

5.2. Source south-west of [PMH2004] 3 in the M 33 
field 

This object, named 3b in Table 1, was not included in 
the catalogue of X-ray sources of PMH2004. Its position 
on a EPIC PN CCD boundary, about 30" south-west of 
[PMH2004] 3 and far off-axis, prevented its detection by 
automatic procedures. However, the 0.2 1 keV X-ray im- 
age (see smoothed contour overlay in Fig. 7) clearly indi- 
cates X-ray emission, most likely originating from the M3 
star USNO-Bl 1204-0019500 (R magnitude 14.66). 

The spectrum we obtained for this object is that of an 
M star. 
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Table 2. X-ray properties of the objects of Table 1. 
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M 31 



ID^ 


U (mW/m^) 


log{fx/fopt) 


HRl 


HR2 


26 


8.3e-15±9e-16 


-3.6 ±0.1 


0.47±0.08 


-0.52±0.09 


31 


2.4e-14±le-15 


-4.0 ±0.1 


0.71±0.04 


-0.60±0.04 


49 


6.9e-15±7e-16 


-4.2 ±0.1 


0.74±0.08 


-0.48±0.07 


59 


7.0o-15±lc-15 


-4.6 ±0.2 


0.39±0.16 


-0.39±0.19 


101 


6.6c-15±5c-16 


-4.1 ±0.1 


0.75±0.05 


-0.69±0.06 


137 


l.lG-15±4e-16 


-4.8 ±0.4 


0.16±0.21 


-0.47±0.27 


168 


2.9e-14±le-15 


-3.6 ±0.1 


0.68±0.03 


-0.58±0.04 


217 


2.0e-15±4e-16 


-5.2 ±0.2 


0.50±0.12 


-0.93±0.13 


464 


1.5c-15±5c-16 


-5.1 ±0.3 


0.24±0.24 


0.07±0.22 


479 


4.9e-15±5e-16 


-4.6 ±0.1 


0.68±0.08 


-0.53±0.08 


498 


2.9e-15±6e-16 


-4.7 ±0.2 


0.72±0.17 


-0.60±0.16 


553 


1.5e-14±le-15 


-4.0 ±0.1 


0.58±0.07 


-0.63±0.08 


615 


2.1e-15±6e-16 


-4.4 ±0.3 


0.30±0.18 


-0.79±0.22 


733 


2.0e-15±5e-16 


-4.5 ±0.3 


0.19±0.23 


-0.33±0.27 


M 33 


3 


6.33c-15±1.7c-15 


-4.1 ±0.3 


0.37±0.21 


-0.64±0.26 


28 


2.58c-14±2.1c-15 


-3.0 ±0.1 


0.44±0.07 


-0.33±0.07 


39 


6.29e-15±8.7o-16 


-2.6 ±0.2 


0.49±0.10 


-0.24±0.11 


54 


7.51e-15±1.4e-15 


-3.9 ±0.2 


0.64±0.15 


-0.75±0.14 


77 


2.34e-14±3.3e-15 


-3.6 ±0.2 


0.67±0.11 


-0.48±0.12 


92 


2.20c-15±8.2o-16 


-5.2 ±0.4 


0.85±0.14 


-0.75±0.18 


122^ 


4.85e-15±1.7e-15 


-3.0 ±0.4^ 


0.94±0.32 


-0.36±0.29 


142 


2.78e-15±7.8e-16 


-5.5 ±0.3 


0.73±0.20 


-0.51±0.22 


182 


2.94e-15±5.4e-16 


-5.9 ±0.2 


0.35±0.13 


-0.54±0.16 


196^ 


1.49e-14±8.6e-16 




0.38±0.05 


-0.65±0.05 


200 


1.43c-14±9.4o-16 


-3.9 ±0.1 


0.61±0.05 


-0.52±0.05 


204 


5.17e-15±6.5e-16 


-4.6 ±0.1 


0.46±0.10 


-0.73±0.11 


206 


6.05e-15±l.le-15 


-3.5 ±0.2 


0.61±0.32 


0.32±0.20 


281 


2.10c-15±5.4e-16 


-5.0 ±0.3 


0.31±0.17 


-0.69±0.18 


297 


7.80c-14±2.2e-15 


-2.5 ±0.1 


0.56±0.03 


-0.10±0.03 


337 


2.09e-14±1.9e-15 


-3.3 ±0.1 


0.33±0.12 


0.30±0.09 


358"* 


9.51e-15±1.6e-15 




0.56±0.20 


-0.14±0.16 


372 


2.60e-15±7.7e-16 


-4.7 ±0.3 


0.80±0.18 


-0.54±0.21 


406 


1.47e-13±5.3e-15 


-3.6 ±0.1 


0.47±0.04 


-0.41±0.03 



^X-ray source, as in Table 1 

^ For the calculation of log{fx/ fopt) for this object wc have used the R magnitude of the USNO-Bl counterpart given in Table 
1, which is most probably a blend of two fainter stars (see Section 5.3). Therefore, the derived ratio log{fx/ fopt) should be 
treated as a lower limit. 

^ This source was not identified in the USNO-Bl catalogue (see section 5.4), but only in the 2-MASS catalogue. Thus, there is 
no estimate of fopt, and therefore of log{fx/ fopt) for this object. 

''The published optical and X-ray fluxes do not necessarily correspond to the same contour of the galaxy, thus we did not 
attempt to calculate log{fx/ fopt) for this object. 



5.3. [PMH2004] 122, in the M 33 field 

This source correlates with the bright USNO-Bl star 1210- 
0020531. However, inspection of the deeper 4m Mosaic 
images of M 33 from the Kitt Peak National Observatory 
archive (from the Local Group Survey of Masscy ct al. 
2002) shows that there are actually two relatively bright 
stars separated by ^6" in NNE/SSW direction which are 
probably blended into one star in the USNO-Bl catalogue. 
We have obtained spectra for both stars. They both have 



G-type spectra, i.e. both stars are compatible with the 
known X-ray properties of the X-ray source [PMH2004] 
122 (see Table 2). Thus, it has not been possible to confirm 
which of the two stars is the optical counterpart of the 
X-ray source. More accurate position of the X-ray source 
would help resolve this question. Unfortunately, the source 
is outside the FOV of the archival Chandra observations. 



8 



Please give a shorter version with: \authorruiming and/or Xtitileriuming prior to \maketitle 



Fig. 8. Optical spectrum (not flux normalized) of the op- 
tical counterpart of [PMH2004] 196 in M 33, identifying 
it as a late M star with H„ emission. 

Fig. 9. Chandra ACIS S hghtcurve of [PMH2004] 196 in 
M 33 during observation 786 integrated over 1000 s. Time 
zero corresponds to the start of the observation. 

5.4. [PMH2004] 196, in the M 33 field 
(2MASS01334186+3038491 ) 

Due to the position of [PMH2004] 196 close to the cen- 
ter of M 33, its optical counterpart was not catalogued 
in USNO-Bl but only in 2MASS. We therefore give the 
2MASS position in Table 1. The optical spectrum identi- 
fies it as a late M type star, which shows significant H„ 
emission (see Figure 8). The emission line is clearly seen in 
all four separate exposures obtained for this star (the com- 
bined spectrum from all four exposures is shown in Fig. 
8). The equivalent width of the emission line is estimated 
to be 5.7 ±0.3 A. The presence of Balmer emission lines in 
late M dwarfs is quite frequent and is linked to variability 
(e.g. Jaschek & Jaschek 1987). Flare stars are expected to 
exhibit a late spectral type (typically dwarf M-type) and 
to present emission in Call (not included in our spectral 
range) and in the Balmer lines, during the quiescent phase. 
The strength of the Balmer emission reaches a maximum 
during the second phase of a flaring event (i.e. after the 
first strong continuum phase). 

The object displays long-term X-ray variability, chang- 
ing its average flux by a factor of 6, most likely due to dif- 
ferent amounts of flaring activity during the observations. 
A long term XMM- Newton EPIC light curve is presented 
in Misanovic et al. (2005, their source 170). The source is 
also found to be variable in the analysis of Chandra data 
(Grimm ct al. 2005). We analyzed Chandra data from the 
archive to search for flaring. While the source was faint 
during observation 1730, it showed strong flaring activity 
during observation 786, when it was found to be brighter 
on average (see Figure 9). Although the source is faint in 
X-rays we managed to produce an XMM-Newton EPIC 
X-ray spectrum, by combining the counts from three ob- 
servations. We have obtained an acceptable fit to the spec- 
trum for a plasma with kT of 2.8 keV. Such a high tem- 
perature is expected during flares. 

5.5. IPMH2004] 206 in the M 33 field (USNO-Bl 
1209-0020969) 

This is an early-K star which has a somewhat high hard- 
ness ratio and it was classified as a possibly hard source by 
PMH2004. The X-ray flux does not show signiflcant vari- 
ability within the errors, which are quite high due to the 
faintness of the source. Moreover, the optical spectrum 
of the optical counterpart shows no indication of emis- 
sion lines that might indicate flaring activity. However, 



Fig. 10. Hardness ratio (HR2 as defined in PMH2004) 
as a function of epoch of observation for [PMH2004] 206, 
which corresponds to source 176 in Misanovic et al. (2005). 
The time is given in days, with the reference day at JD 
2450814.5. 

the hardness ratio of the source does show some evidence 
of variability between observations that are spread over 
three years (Figure 10). As can be seen in this figure, the 
source often appears to have a much softer hardness ratio 
than that given in Table 2 (and plotted on Fig. 6). 

There seems to be at least one other early K giant that 
shows significant coronal activity and some temperature 
variability (Audard ct al. 2004). The authors suggested 
that most of the coronal heating in that star was con- 
tributed by small flares. 

5.6. [PMH2004] 297, in the M 33 field (USNO-Bl 

1209-0021478) 

In PMH2004 this object was misidentifled as a back- 
ground elliptical galaxy (object MD 53 in Christian and 
Schommer 1982). It is actually a stellar source to the 
North of this galaxy. Indeed, the spectrum is a stellar 
spectrum of late G-type. 

The source displays significant long-term variability 
in X-rays, changing its flux by a factor of 4. The light 
curve is presented in Misanovic et al. (2005, their source 
253). [PMH2004] 297 is one of the two sources among the 
objects studied here for which an X-ray spectrum could 
be produced. Preliminary analysis of the spectrum indi- 
cates that it is typical of very active stars (for example 
see Guedel 2001, 2004). We have obtained an acceptable 
fit with a combination of two temperature plasmas with 
kT of 0.8 and 2.3 keV. The high temperature component 
is most probably produced by flaring or some other type 
of stellar activity, which may also explain the significant 
X-ray variability of this source. 

5.7 [PMH2004] 337, in the M 33 field (USNO-Bl 
1204-0021198) 

This object has a hardness ratio (HR2)^ value that is too 
high compared to other stars in the sample with sim- 
ilar X-ray to optical flux ratios. The optical spectrum 
of the star shows no indication of emission lines that 
might indicate flaring activity. However, the X-ray source 
itself is relatively close (0.99 arcmin) to a very bright 
hard source, [PMH2004] 335 (HR2= 0.19 ± 0.01), which 
may disturb the count rate and hardness ratio determi- 
nation, as it is two orders of magnitude more luminous 
(/^ = 1.56e - 12 ± 1.6e - 14) than [PMH2004] 337. 

^ The source is covered by three individual XMM-Newton 
observations and detected above the selected threshold in two 
of them. The HR2 seems to be approximately the same (0.24 
and 0.28) in both observations separated by several months. 
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Fig. 11. Images of the galaxy LED A 5899, in U (top 
panel) and V (bottom panel). The circles show the X- 
ray positions of the source and the corresponding 3 sigma 
error boxes from PMH2004 (bright circle) and Misanovic 
et al. (2005 dark). 

Fig. 12. Optical spectrum (not flux normalised) of Galaxy 
LEDA 5899. 

5.8. [PMH2004] 358, in the M 33 field (LEDA 5899) 

This source corresponds to the spiral galaxy LEDA 5899 
(CGCG502-114) which is found within the M 33 D25 el- 
lipse and shows a somewhat disturbed appearance. Figure 
11 shows images of the galaxy in U and V (top and bot- 
tom panel), extracted from the 4m Mosaic images of M 33 
from the Kitt Peak National Observatory archive (from 
the Local Group Survey of Massey et al. 2002) . The over- 
layed circles show the X-ray positions of the source and the 
corresponding 3 sigma error boxes from PMH2004 (red cir- 
cle) and Misanovic et al. (2005 yellow circle). The optical 
spectrum of the galaxy that we obtained here is presented 
in Figure 12. Some characteristic lines are marked. Our 
determined redshift (0.034) agrees within the errors with 
the value of 0.035 given in Huchra et al. (1999). X-ray 
variability could not be established for this object, due to 
lack of sufficient data. The X-ray luminosity of the galaxy 
in the band (0.2-4.5)keV is found to be 2 - 3 x W^^eig 
s~^, a luminosity typical of normal spiral galaxies (see 
e.g. Fabbiano et al. 1992). 

6. Summary 

We have obtained optical spectra for the 14 brightest op- 
tical counterparts of the XMM-Newton source catalogue 
of the M 31 field (PFH2005) and the 21 brightest counter- 
parts of the XMM-Newton source catalogue of the M 33 
field (PMH2004), using the 1.3m Skinakas telescope in 
Crete, Greece. All of the M 31 sources and all but one of 
the M 33 sources are foreground stars, of spectral types 
between A5 and M5. One of the M 33 sources (lying within 
the D25 ellipse) is a background galaxy. The majority of 
the objects have X-ray hardness ratios and X-ray to opti- 
cal flux ratios (see Table 2) that are consistent with being 
a foreground star. One of the stars close to the M 33 cen- 
ter is a late M dwarf with Hq emission, probably a flare 
star, also displaying strong X-ray variability. 
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